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SUMMARY 



This report presents the results of a combined theo- 
retical and experimental investigation conducted at the 
Georgia School of Technology on the static thrust of the 
lifting airscrew of the type used in modern autogiros and 
helicopters . 

The theoretical part of this study is "based on 
Glauert ! s analysis hut certain modifications are made that 
further clarify and simplify the problem. Of these changes 
the elimination of the solidity as an independent parameter 
is the most important. 

The experimental data were obtained from tests on four 
rotor models of two, three, four, and five blades and, in 
general, agree quite well with the theoretical calcula- 
tions. 

The theory indicates a method of evaluating scale ef- 
fects on lifting airscrews, and these corrections have 
been applied to the model results to derive general full- 
scale static thrust, torque, and f igure-of-mer i t curves 
for constant-chord, constant-incidence rotors. 

Convenient charts are included that enable hovering- 
flight performance to be calculated rapidly. 



INTRODUCTION 



The problem of greater safety in flight is today com- 
manding more and more attention. Two different methods of 
attack are being developed at present. One of these con- 
sists of improving the conventional fixed-wing airplane 
through such modifications as Handley Page slots, wing 
profiles giving smooth maximum lift characteristics, meth- 
ods of obtaining more complete rolling and yawing control 
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in stalled flight, and other special devices. The alterna- 
tive method is that of developing a type of aircraft in 
which there will always oc relative motion "between the lift' 
ing surfaces and the air, regardless of the motion or atti- 
tude of the aircraft as a whole. This type is exemplified 
by the autogiro and the various experimental helicopters, 
of which the Br egue t -Do rand is the most outstanding recent 
example (reference l). 

In order to investigate the possibilities of the 
r o t at ing-wiag type of aircraft, a general study of the ver- 
tical notion of the lifting airscrew has "been undertaken 
at the Daniel G-uggenhe im School of Aeronautics of the 
Georgia, School cf Technology. This project is receiving 
f inanci 1 smp#$rt from the National Advisory Committee for 
Aeronautics and the State Engineering Experiment Station 
of C-oor g i a . 

The purpose of this re;oort is to present the results 
of the first part of this investigation, which covers the 
phase of static thrust or hovering flight of the helicop- 
ter. C-iauert ! s assumptions (reference 2) furnish the 
"background for the theoretical portion of the study. Hew- 
ever, the induced velocity through the rotor is determined 
on the oasis of vortex theory rather than "by using the con- 
cept of the "actuator disk." This change has "been made 
hecause the vortex theory offers a much clearer picture of 
the mechanism of airscrew thrust without materially com- 
plicating the derivation of the induced velocity equation 
which is identical for both methods. 

The experimental part of the analysis provides numer- 
ical values of such parameters as are essentially empirical 
and serves to show the agreement "between the calculated and 
actual values of thrust and torque for four different ro- 
tor mo lei s . 

STATIC THRUST ANALYSIS 



Ba^ie __a s r.umptions In treating the complex prohlem 
of the lifting airscrew, it is necessary to make the fol- 
lowing simplifying assumptions : 

1. The number of "blades may "be taken as infinite. 

2, Induced angles of attack are sufficiently small 

so that the angle may he suhstituted for the 
sine and tangent, and the cosine replaced by 
unity . 
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3. Rotational and radial components of velocity and 

tip effects may "be neglected. 

4, The slipstream contraction may "be neglected. 

Induced velocity*** The first step in the static thrust 
analysis is to determine the downward velocity induced hy 
the motion of the "blades. In the simplest case, each "blade 
of the rotor may "be replaced "by a rotating lifting line from 
the ends of which spring trailing vortices. The tip vortex 
will obviously form a helix in space while the vortex ema- 
nating from the blade root will he concentric with the axis 
of rotation, and in the i^resent discussion its effect will 
he neglected. 

Since we are assuming an infinite number of blades the 
problem consists of determining the velocity induced nor- 
mal to the plane of the rotor by a cylindrical surface of 
vorticity, bounded on one end by the rotor and extending 
downward to infinity. Figure 1 depicts this cylinder with 
its top in the xy plane and its axis coincident with the 
-2 axis. The ring s of width dz is an element of this 
surface normal to the axis. 

If we take the total circulation produced by the rotor 
blades as F t the circulation strength of the ring s 

d-T 

will bo ™ dz, and by neglecting higher order terms the 

CL Z 

ring becomes equivalent to a circular vortex element of 

radius R. It will be noted that §*- = constant. 

dz 

Now the potential at a point P due to a closed vor- 
tex element may be expressed in the following form (refer- 
ence 3 ) : 

where r is the circulation about the element 

and uu is the solid angle at P subtended by the 

element 

Thus, if the point P be in the plane of the rotor, 
as shown in figure 1, we may write: 
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where im is the solid angle at P subtended by, the vortex 
ring s. The total potential at P will then "be 



dz 
4tt 



<t' dz 



(1) 



The velocity at P, due to s, normal to the rotor 
plane, may now "be obtained "by differentiating with, respect 
to z : 

- r r 

-p = |f - -ft- U(z 2 ) - ■»(«*> 



and for the entire cylinder where 



uu(Z 2 ) = 0 and w(Z l ) 



~ oo and Z 1 



2tt 



0 : 



i fir 

. = - « = constant 
P 2 dz 



(2) 



It is important to note also that at the point p 1 in the 
plane of the rotor "but "beyond the blade tips, uu(Z 2 ) = 

tt)(Z x ) = 0 and 



0 



(3) 



Thus v;e have established the fact that the vertical com- 
ponent of induced velocity is constant over the rotor 
disk since V = constant. Moreover, outside the disk this 
component is zero. It should also "be observed that, if 
the point ? is moved downward inside the vortex cylinder 
to a groat distance from the rotor disk, the solid angle 
subtending it will become 4tt , since the cylinder may 
then bo considered doublv infinite in extent. Consequent— 



"D 



= - d T 

dz 



(4) 



thus corroborating- the well-known fact that the slipstream 



velocity in the plane of an 
large di stance downstream. 



airscrew is half that at a 



To investigate the general case of the induced veloc- 
ity due to any distribution whatsoever of the circulation 
along the blade, we may proceed as follows. Consider a 
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"blade element of length dr at a distance r from the 

axis, as shown in figure 2. Let the circulation at r "be 

I\ The circulation at r + dr will then be (T + |£ dr). 

dr 

According to equation (3) the induced velocities out- 
fii^e the respective cylindrical surfaces will be zero, 
whereas An_s_ide they will be 



d (r + |r dr ) 

w = - i &£ 

2 dz 



- — from 
2 dz 



to r + dr 



and 



1 

2 [dz 



dz 



= 0 



from 



these results being obtained by neglecting the higher order 
infinitesimal s . 

This simple demonstration verifies the independence 
of airscrew "blade elements, which is a customary assumption 
in modern airscrew analysis. 

Having obtained the induced- vel oci ty relationships due 
to the vortex field, we will now derive a general expres- 
sion for the circulation F as a function of the blade- 
profile characteristics. Consider an element of a blade of 
length dr at a distance r from the axis of rotation. 
The thrust on this element will "be 



dT = p - 0 r dr = - C 7 

B 2 1 



Q 2 r 2 c dr 



(5) 



where 



B is the number of blades 

c , blade cho r d 

fi, angular velocity 



Making the customary assumption that varies linearly 



with angle of attack 
ure 3 : 



a , 



we may write with the aid of fig- 



a o a 



a Q (6 - cp ) 
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where a Q is the lift curve slope for two-dimensional 
f 1 ow 

6, blade incidence angle 
i w w 

cp , tan = — — = incuced angle of flow 

\ L I \ L -L 

From equation (5) re may now write the circulation as 

r = ~5 ao ( 6 - cp) Q r (6) 

The distance "between successive turns of the vortical he- 
lix is 



2l[r 
B 



sin cp 



2rrr 
B Qr 

- 2ttw 
n 3 ' 

and. since — = — — = 2w 

d z B.d 

we ootain from (6) that 



a 0 c 



(6 - cp) fir 



2w 



2ttw 

n 3 



a„ Be 



or w = - 



o _ (G - cp) n 2 r (7) 



8 TT 

In this analysis we shall confine ourselves to the case of 
hlad.es of constant chord and. hence the rotor solidity, as 
usually expressed, is 

_ Be R Be 
~ TT E 2 T ttR 

Putting x m | 
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equation (7) "becomes: 



o r 




The useful root of this quadratic is 



qp 




Dividing "both sides "by o* we obtain: 




and putting 



qp 



cp 
cr 



we have 

: % = 2* L ✓ VTV ' * W a m " 8 j 

which is. identical with the result obtainable "by the use 
of the momentum equation, i.e., the "actuator disk" method* 

It should "be noted that equation (8) is a general ex~ 
pression for the induced flow angle in terms of the radius 
variable x and the two parameters 6 a and a 0 . However, 

a 0 is substantially constant, varying but slightly with 
profile thickness and Reynolds Number. Thus, equation (8) 
effectively gives cpQ- as a function of the single parame- 
ter 6 a , which results in a marked simplification of the 
analysis as originally presented by G-lauert (reference 2 ) • 

In the first part of the previous analysis each blade 



/ 



a, 



i A 



/a 



a. 



( a \ 



8 



N.A.C.A. Technical Note No. 626 



was replaced by a lifting line, i.e., F = constant 
Therefore, equation (6) may he written: 



r = ~£ a Q (6 - cp) Qr 



Be 



2 -o ( 0 o - <P 0 > QR 



or 



(6 - cp) r = (G 0 - cp n ) R 



where 8 Q and cp Q refer to the tip of the "blade, hut 
since cp = w/fir 



e = 



X 



(9) 



This is the variation in G along the "blade that is re- 
quired to give constant circulation and uniform induced 
velocity. Over the effective portion of the "blade 

10 6 
x > - > £ = as tan"" (anDrox.) 
4 x x 



which defines the screw surface described "by each hlade, 
if it were moving in a solid medium. Therefore, equation 
(9) may he considered to represent the "constant pitch 11 
airscrew, and the induced angle function becomes 



+ 4 - 



9^ 2. 



■ 7 



(10) 



For the type of rotor with untwisted "blades, i.e., constant 
inc idence 



'a 



x - - s 



(11) 



where 0- = constant 

Thrust The thrust produced hy the rotor may he ex- 
pressed as follows: 

T = E tt R 2 Q 2 R 2 Cm 

2 1 
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where C T is the nondimensional thrust coefficient. 

To determine the thrust in terras of the rotor parame- 
ters, we write : 

1 

T = | Be j C L Q 2 r 2 dr 



P 
2 



3c R 3 a n / (6 - cp) n 2 x 2 dx 



o r 



C_ = a a / (G - cp) x 2 dx 



and dividing "by cr c 

C m 



(6 a - cp a ) x 2 dx 



(12) 



This new thrust coefficient will he designated as 



For the constant Ditch rotor 



where qp = \ [ J (*f) * 4 (-/) 6^ - 



(13) 



And for the constant-incidence rotor 



I as a. 



_1_ 

2x 



x 2 dx 
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and upon integrating 



j _ a /?£ + I _ 1 Al _ £ + 

^ ~~ 0 \ 3 4 V g / „ 2 5 V 8 > 3 



5 / -.3 



" ^ (¥) 



(14) 

15 V 8 / j J 



Torq ue." The rotor torque may "be written as 
Q = | tt H 2 Q 3 R 3 C Q 

Cq = nondimensional torque coefficient 

The torque may "be divided into three parts analogous to the 
partition of drag on an airfoil. These are: 

1. Induced torque. 

2. Minimum profile torque. 

3. Profile variation torque. 

The induced torque due to the inclination of the air- 
foil lift vector, figure 3, is 

E 



O ui = | 3c / C L sin cp n 2 r 3 dr 



or 



= g Be R 4 a 0 / cp (6 - cp) x 3 dx (approx.) 
o 

C 0 = 0" a 0 / cp (6 - cp) x 3 dx 



0 



The minimum profile torque is 



N.A.C.A. Technical Note No. 626 



11 



Q,c = | le / 5 r 3 dr 



2 



o 



l 

P Be R 4 6 / ft 2 x 3 dx 



2 



./ 

o 



or C n = r 6 / x 3 dx 



where 6 = 

o . 

mln 

Assuming that the profile drag variation may "be ex- 
pressed as 

A C-n = € a 2 

o 

where e a const. 

The profile variation torque may "be written as 



Q t = £ Be / e a 2 Q 2 r 3 dr 



'J 2 



l 



| Be R 4 c / (6 - cp) 2 n 2 x 3 dx 



• / 

0 



1 



or 



P 2 

Cq = a c y (e - cp) x 3 dx 



°0 



/ 

0 

?or the constant pitch rotor the total torque coeffi- 
cient may he expressed in the new form to eliminate C as 
an independent parameter as follows: 



- A l ,10 2 

Q a . Jt / x3 dx + € /" - cp a ) x* dx * 

° ,^ e 0 

+ a 0 / cp r (-^ - cp r ) X 3 dx (15) 

0 



12 N.A.C.A. Technical Note No. 626 



But since 




In like manner the torque coefficient may "be obtained- 
for the constant incidence rotor as 




(17) 



where, as in the thrust equation (14) 




Figure of merit Since for the condition of static 
thrust the rotor is not moving in translation, the ordi- 
nary concept of efficiency must "be modified in order to de- 
vise a method of measuring the "lifting efficiency" of the 
helicopter airscrew. G-lauert in reference 2 has designated 
this criterion as the "figure of merit" and has defined it 
a s 



n 3 /2 
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This expression may "be arrived at "by noting that 

2 2 

iji C rj\ Cl R 

P = C Q Q 3 R 3 

where P is the power applied to rotor; and to eliminate 
the tangential velocity QR, we may write: 

t-'JL °i 3/ * 
____ . __ 

and in terms of the new coefficients 

U> . . -°— (18) 



0 0 \ Qct 



^a 5 > 

The theoretical maximum or ideal value of M would 
occur if 8 and e were zero and its value may he ob- 
tained easily for the constant pitch rotor. Referring to 
(7a), 



T CT = - ? V ' w ° -^V = 4cpJ 

2 °o 



Q _ 21 = 



o r 



m 3 / 2 

1 a 



and hen ce M T - = 2 

The ideal figure of merit for the constant incidence 
rotor is not a constant and cannot "be easily obtained as 
an analytical expression. However, the numerical values 
have "been computed and these are given later. 
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The figure of merit may be defined on the "basis of uni- 
ty, as in the usual definition of efficiency by writing 



m 3/2 
1 GT 



M 



and this form will "be used in the subsequent discussion, 



EXPERIMENTAL LATA 



In order to obtain the values of the a.irfoil parame- 
ters appearing in the theory and to enable the calculated 
valiies of the rotor coefficients to be compared with tests 
on actual rotors, two sets of tests were made in the 9- 
foot wind tunnel of the Daniel Guggenheim School of Aero- 
naut teg* 

Rotor model tests.- Four rotor models having two, 
three, four, and five blades and a diameter of 5 feet, 
were tested. Five blades and three hubs were used. The 
blades were identical and interchangeable, thus making 
possible the four rotor combinations. 

The blades had the IT.A.C.A. 0015 symmetrical profile 
(reference 4") and were of 2-inch chord from the tip to a 
radius of 5 inches. From this point they were faired into 
a ?/4-inch circle at a radius of 1.5 inches. The diameter 
of the hubs was 3 inches. Horizontal hinges were located 
in the hub at a radius of 1 inch to permit vertical artic- 
ulation of the blades. A split sleeve with clamping screws 
at the inner end of each blade enabled the blade incidence 
to be changed. The blade plan form was straight and there 
was no twist. 

In order to prevent any tendency to twist when in op- 
eration, the "blades were statically balanced about the 
quarter-chord point by using a steel leading edge recessed, 
into the laminated wood of the blade and the quarter-chord 
line was radial. The blade tip in cross section was semi — 
c i rcular • 

The various rotors were mounted on the frame of the 
wind-tunnel balance as shown in figure 4. The rotor drive 
shaft was horizontal as shown, and was 12 inches long, 
terminating in a small gear box*'at the center of the tun- 
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nel jet. This gear hex was held hy a torque tube support- 
ed on a tripod attached to the "balance frame and extending 
below the jet to a motor-driven worm and pinion, enabling^ 
the angle of attack of the rotor to he varied through 360°. 

Inside this torque tuhe was the main drive shaft which 
was turned hy a 1-hor sepower 3-phase induction motor 
through a V-belt drive. Since all this equipment was 
mounted on the six- component "balance frame, it constitued 
a dynamometer capable of measuring the thrust and torque on 
the rotor for any position of the model. 

The six-component balance of this wind tunnel upon 
which both the rotor and airfoil tests were made, is some- 
what unusual in that it utilizes a hydraulic method of 
transmitting the forces to a central point where they are 
automatically balanced and indicated pneumatically. The. 
U frame which holds the model is held in place by six tubu- 
lar struts, each of which terminates in a hydraulic cell. 
These cells consist of a shallow cylinder containing a 
loose-fitting piston. The l/8-inch annular space between 
the piston and cylinder is covered with 0.020-inch thick 
pure gum rubber, thus providing a tight but substantially 
f rict ionless seal. The cells are completely filled with 
distilled water and the pressure due to the force applied 
on each cell is transmitted through small copper tubing to 
identical cells mounted on the frame shown in the right 
background of figure 4. Here the forces are combined by a 
rigid connection between the pistons of the appropriate^ 
cells so as to give the total force. These forces are in 
turn applied to the pistons of similar cells actuated by 
air, and the pressures in these air cells are controlled^ 
automatically hy means of special piston-type valves^which 
apply either pressure or vacuum to the cells, depending 
upon the direction of the forces. Since the pressure in 
each air cell at any time is a measure of the particular 
force applied to it, these cells are tapped and the pres- 
sure transmitted through rubber tubing to a manifold which 
is fitted with six valves. From this manifold another rub- 
ber tube goes to an additional cell unit of the same type 
mounted on the platform of a dial scale. Thus a single 
operator may read all the forces in any desired sequence 
by opening the appropriate manifold valves. The sensitiv- 
ity of the balance can be changed by substituting balanc- 
ing air cells of different sizes or by using water or al- 
cohol manometers in place of the dial scale. Forces as 
small as 1 gran or as large as 150 kilograms can be meas- 
ured with the present arrangement. This range could easi- 
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If "be extended. The system is practically null in opera- 
tion, since the maximum movement of the balance frame is 
never greater than about l/64 inch. 

The balance has proved quite satisfactory over a peri- 
od of three years, the only attention necessary "being an oc- 
casional replacement of the rubber portions of the dia- 
phragms and cleaning of the piston valves, It has proved 
to be convenient and flexible in its operation, since it 
not only enables a single operator to make all ordinary 
tests but it is also possible to move the indicating mech- 
anism wherever desired. 

The average rotational speed of the model rotors was 
approximately 960 r.p.m. This varied slightly with the 
number of blades and with the incidence because of the in- 
duction motor slip and the creep of the belt drive under 
load, but not sufficiently to cause an appreciable change 
in Reynolds ITumber. The speed readings were obtained by 
means of an ordinary hand tachometer and stop watch, the 
average of thr«ee 30-second readings being used. 

In order to obtain consistent and reliable thrust and 
torque readings, it was found that certain precautions had 
to be taken. It was necessary to set the blade- incidence 
angles within plus or minus 0.05°. This was finally accom- 
plished with the aid of the apparatus shown in figure 5, 
The rotor hub was clamped rigidly, to a horizontal support 
and the tip of each .blade in turn was attached by a special 
clamping device to a sensitive inclinometer mounted on ball 
bearings and count erwe ight ed so that it imposed no restraint 
on the flexible blade. The inclinometer was sensitive to 
within one minute of angle. 

There was found to be a slight lack of uniformity in 
the blades due to warping. Since, for uniform results, it 
was necessary to set the blade angle with respect to the 
zero thrust angle, the thrust curves of each blade were ob- 
tained by running them singly with a "balancing counter- 
weight. Two of the blades showed no effective twist and 
the worst one had a twist of 40 minutes. In making the fi- 
nal tests these twist corrections were used in setting the 
blade angles. 

Various positions of the roodel with respect to the 
tunnel were investigated and the final position chosen was 
that shown in figure 4. with the slipstream toward the left 
since in this position no interference effects were appar- 
ent. 
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The precision of the measurements nas as follows: 



Thrust 



±1 percent 



Torque 



±1 percent 



Minimum torque 



±2 Dercent 



Rate of rotation 



Blade incidence angle 



±J percent 
±0.05° 



Airfoil tests *- In the above mathematical analysis, 
the airfoil profile characteristics of the rotor blades 
were assumed to "be as follows: 



where a is the effective angle of attack of the airfoil. 

In order to obtain quantitative results from the the- 
oretical equations of thrust and torque, it was necessary 
to determine specific values of a and Oj) , "both of 

°min 

which are functions of the Reynolds Number. This was done 
by testing an airfoil of N.A.C.A. 0015 profile in the wind 
tunnel at approximately the same Reynolds Number (242,000) 
as that of the rotor-blade tips. This airfoil had a span 
of 6 feet and a chord of 6 inches and was mounted in the 
wind tunnel as shown in figure 6, which is a view of the 
set-up looking in the upstream direction. 

The model was supported on a streamlined fork located 
at the midspan and quarter-chord point and on a small oval 
rod which was connected to a short sting attached to the 
trailing edge. This latter rod was actuated by a push- 
pull rod sliding inside the torque tube which held the 
fork and was attached to the balance frame. The push rod 
was actuated by a second motor-driven worm and pinion mech- 
anism permitting remote control of the angle of attack* 
The angle readings were effected by means of a simple di- 
rect-current bridge system with the galvanometer and the 
adjustable leg of the bridge mounted on the balance -control 
table • 






+ € a 2 
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The "balance frame and the model- supporting tripod were 
shielded as shown to reduce the tare forces. These tare 
forcer, were measured with the model supported at the tips 
and with the fork and incidence rod in ^lace "but not touch- 
ing the airfoil and were found to he zero for the lift and 
about 60 percent of the minimum profile drag. 

The air flow in the jet is reasonably smooth and uni- 
form owing to the 5 to 1 contraction in the collector. 
The variation of dynamic pressure over the span of the mod- 
el is within plus or minus 1 percent. A3 previously men- 
tioned, the jet is 9 feet in diameter and 12 feet long, » 
and the model is located with its quar tor-chord point on 
the lateral axis, of the "balance, which is 3 feet downstream 
from the entrance cone. 



In these tests the dynamic pressure was held constant 
by means of an alcohol manometer connected to four static 
plates located in the large section of the entrance cone, 
the static plate pressure having been previously calibrat- 
ed against pi tot surveys nra&e in the model position. 

The precision of measurements in the airfoil tests 
was as follows : 



Lift .... 
Drag . 
Minimum drag 
Velocity . . 



±1 percent 
±1 percent 
±2 percent 
±h Percent 



Anrle of attack ±0.05° 

Reduction of data*- The formulas used for reducing J :he 
thrust and torque measurements from the four rotor tests 
are as follow?. : 



2 _2 2 
O TT R Q R 



-a 



V 



CF 



'P TT H 



8£\ 
4 /' 



' W 1 4 a 2 / 
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The thrust and torque coefficients and Cq 1 , re- 

spectively, are plotted against angle of "blade incidence 
in figures 7 and 8. 

Reduction of the data from the airfoil tests requires 
an interpretation of the wind-tunnel jet "boundary effects, 
A study of the data revealed that the interference of the 
"balance shielding was appreciable and, in an attempt to 
evaluate this interference, the work of Tani and Taima 
(reference 5) on the "boundary influence of partial enclo- 
sures consisting of circular arcs, was consulted. It ap- 
peared that the "balance shielding had an effect equivalent 
to an arc enclosure of approximately 145°, which results 
in o. zero correction for induced drag and angle of attack 
"but a balance alignment correction of about 0*30° due to 
the upward inclination of the air stream produced by the 
boundary. This alignment correction was, therefore, ap- 
plied to the original force test data. The results were 
then corrected to infinite aspect ratio by the customary 
formulas (reference 6). However, a small asymmetry re- 
mained, the inverted tests giving slightly higher values 
o.f lift and drag, and it was therefore necessary to draw 
an average curve between the points. Figure 9 shows 
these final curves with the minimum drag coefficient sub- 
tracted from the total profile drag to give the. profile- 
drag variation curve required for the determination of € # 

COMPARISON OF THEORY AND EXPERIMENT 



Thrust.- The experimental values of the new thrust 
coefficient T a are plotted against 6 a in figure 10. 

This figure shows at once that the assumption of an infi- 
nite number of blades and the elimination of a as an ^ 
independent parameter are fully justified since the points 
for the four rotors of solidities 0.0424, 0.0636, 0.0849, 
and 0.1031 all fall very closely on a single curve with 
the exception of the high incidence values where the pro- 
file lift curve no longer approximates a straight line. 

In figure 9 it will be seen that a mean value of the 
lift-curve slope for infinite aspect ratio is 5.75 per 
radian. Using this value in equation (14) the theoretical 
thrust curve, also shown in figure 10, was obtained and 
indicates a fairly good agreement between theory. and ex- 
periment except for small values of Bq* 
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In order to obtain a single curve to define 
the static torque characteristics of the helicopter, it is 
evident from equations (16) and (17) that we must plot 

Q(T ! = lQ,(j - — "§ ) vs. 8<j, In figure 11, therefore, the ex- 
v 4 CT ' 

perimental values Qq- 1 are plotted against 6 a and here 

again we find that the points for the different rotors 
fall on a single curve with the exceptions noted above for 
thrust • 

In order to obtain the theoretical curve, it is nec- 
essary to determine the value of € , the coefficient of 
profile-drag variation. This value, as shown in figure 9, 
is 0.75 at the Reynolds Number of the blade tips. This 
figure shows also that the assumed parabolic variation of 
profile drag is justified, at least for this R'eynolds 
Numbe r . 

Using the value of 0.75 the' torque coefficient curve 
was obtained as shown in figure 11. This curve falls be- 
low the experimental curve, but by increasing e to 1.25, 
excellent agreement is obtained. 

A careful analysis of the assumptions has revealed 

that the only one which night account for the thrust and 

torque discrepancies is that of neglecting the slipstream 

contraction but further study is required on this point. 

As a matter of interest, the experimental values of 
thrust and torque are plotted against each other on loga- 
rithmic paper in figure 12, and indicate that the follow- 
ing simple relation holds quite closely for the model 
tests: 

Q a ' = 0.72 T a 3/2 

The thrust and torque coefficients are also plotted on 
rectangular coordinates in figure 13. These results show 
that the theoretical values of thrust are constantly higher 
than the experimental values although the curves are of the 
same general shape. The approximate equation gives excel- 
lent agreement with the experimental curve up to a Qq- T of 

about 5. Beyond this point the approximate values of 
thrust exceed the experimental values by an increasing 
amount . 
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ANALYSIS OF PARAMETERS 



In the past the interpretation of helicopter airscrew- 
model data has "been made difficult and uncertain "because 
of large scale effects. The magnitude of these differences 
due to scale may "be judged from figure 14, which gives the 
curves of figure of me-rit M vs. "blade-angle factor 6q- . 
for "both model and full-scale conditions. The solidity 
chosen a = 0.06, is approximately a median value for ex** 
isting helicopters and autogiros, and the method used to 
obtain the data for these curves will become apparent in 
the subsequent discussion. It should "be noted that the 
Reynolds Number range of the blade tips of existing heli- 
copters and autogiros is approximately 2.5 x 10 6 to 
3.5 x 10 6 . 

Fortunately, the theoretical portion of this investi- 
gation has revealed the existence of certain parameters 
that are functions of scale and the values of these corre- 
sponding to any Reynolds Number may be obtained from suita- 
ble airfoil tests. These airfoil parameters are: 

a is the slope of lift curve for two-dimensional flow. 

8, minimum profile-drag coefficient. 

€, profile-drag variation coefficient 

Lift-curve slope, a 0 . - For two-dimensional flow the 

slope of the lift curve varies slightly with the Reynolds 
Number and also with thickness and camber as shown in ref- 
erence 4. However, inspection of this reference and fig- 
ure 9 of this report, indicates that a 0 = 5.75 very 
nearly for both model and full-scale Reynolds Number. More- 
over, for different airfoil profiles within the convention- 
al thickness and camber ranges, the variations in a 0 are 
negligible. Thus, no appreciable error may be expected in 
assuming a 0 a const. = 5.75. 

This conclusion enables us to derive the full-scale 
curve of thrust coefficient T a vs. blade- inc i den ce fac- 
tor 6 a . However, it should first be recalled that figure 

10 shows an appreciable discrepancy between the theoreti- 
cal and experimental curves for the constant incidence 
type of rotor for £ < 3 r 5 . Consequently, for the sake 
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of accuracy, 
over 

T 



the faired experimental 



'(J 



this range. For 
- , equation (14 ) must he 
stalling of the model 



> 



3.5 

used 



the 



values should he used 
theoretical values of 



lince "beyond this point 



rotcr "blades 
rotor of solidity a = 0.0424. In 
use of the theoretical values, one 
the single-blade rotor tests, a - 
figure 10, and this falls very 
curve at 6 = 4/93. 

efficient curve thus 
coordinates since 



hegins to occur for the 
order to justify this 
experimental point from 
0.0212, is- plotted in 
falls very close to the theoretical 
In figure 15 the modified thrust co- 

ohtained is plotted on logarithmic 
this form gives approximately constant 



accuracy for reading values from the curves. 

As the incidence of the hlades is increased a point i 
finally reached where the hlades. "begin to stall. In the 
experiments this condition could easily he detected hy the 
loud roaring noise that resulted. However, even "before 
the stall occurs the lift curve departs sufficiently from 
a straight line to cause appreciahle divergence from the 
theory and it is, therefore, necessary to determine this 
limiting value of hla.de incidence. Referring to equation 
(7a), we note that 



Cp" a ( G - Cp) 



8x 



= ( 6 ~ a) 



o r 



^2 a a ^ 

6 - 2E a + a 2 - - = 0 

x 8 



and 



G = a 



+ J ; 



a r 



Inspection of the derivation leading to equation (8) 
indicates that, for hoth constant-pitch and constant-inci- 
dence rotors, stalling will first occur at the hlade tips, 

Hence, the desired limiting value of tip 



i.e., x 

inc i dence 



i s 



a + 



w here 



a 



is the absolute angle of attack (radians). 



A study of the airfoil characteristics given in ref- 
erence 4 shows that the value of a at which the profile- 
lift curve departs noticeably from a straight line, is 
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ahout 14° for practically all airfoils having thickness 
ratio not less than 0.09. Thus, for present full-scale 
lifting airscrews," we may as sun e that 



o r 



0.25 + 

0*25 

cr 



0.25 



0.25 a r 



(7 



The limiting values of 6^ derived from the atove 

o 

equation are sh<*wn in figure 15 for different solidities 
and are also indicated in all the f igure-cf-neri t curves. 

P rofile -drag variation coefficient . € . - In presenting 
the experimental torque data in the form of the new coeffi- 
cient Qq- ! » mention has already "been mo.de of the fact that 



a' = (la - f~ s ) ■ f (*o* 6 a) 



But Qo- r is composed of two parts as indicated in equations 
(16) and (17). These are the "induced torque" and the "pro- 
file-variation torque," and in the latter the quantity € 
appears as a coefficient. The value of this coefficient 
for the rotor-model tests was found to "be 1.25 as shown 1 in 
figure 11, although the airfoil tests indicated a value of 
0.75. (See fig. 9.) This discrepancy is partly due to the 
lower Reynolds Numher at which the inner portions of the 
"blades were operating and probably also to the large inward 
radial velocity components near the tips that have been 
neglected in the theory. To account for this difference, 
we may write: 

c = k c 1 

where e is the profile torque variation coefficient 

t 1 , profile drag variation coefficient 

K = 1*67, correlation factor 

Since no suitable full-scale test data on lifting airscrews 
are known to the writers, we shall assume that K is in- 
dependent of scale. On this oasis the full-scale value of 
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€ may "be estimated from the lower curve of figure 9 which 
was obtained from reference 4 at R.N, = 3.5 x 10 for 
the N.A.C.A, 0015 airfoil used in the rotor models. This 
value of € ! appears to "be about 0.18 and hence e s 0.30 
should he a fair representation of the coefficient for ex- 
isting full-scale rotors. 

In estimating the full-scale torque coefficient Qq- f , 
figure 11 indicates that the theoretical data may "be ex- 
pected to give good accuracy. Hence the theoretical val- 
ues for the constant-incidence rotor' have "been computed 
using c = 0.3 and a Q = 5.75 and yield the torque curve 

Qq- 1 vs. 6<j given in figure 15. 

Variations in c with camber are greater than those 
with thickness "but in general "both are small enough to he 
neglected. 

Minimum profile-drag coefficient, 8 . - This coefficient 

ar>r>ears in the form ~4sr in equations (16) and (17) and is 

4 Q 

obviously constant for a given rotor. The term — is 

the minimum torque and is obtained when the thrust is zero, 
The value of 8 obtained from averaging the four rotor 
model tests was 0.0115 and from the airfoil tests was 
0.0113, a surprisingly good agreement. 

The determination of full-scale values of 8 presents 
some difficulty due to lack of agreement between the re- 
sults of different wind tunnels. Recent large-scale p'ro- 
file-drag tests made at the D.V.L. (reference 7), seem to 
indicate that the values of 8 obtained in the N.A.C.A. 
variable-density tunnel (reference 4) are too high due to 
tip effect and relatively great surface roughness. Until 
such divergences can be reconciled, careful judgment will 
be required in the estimation of 8 for a given airfoil. 

Minimum profile drag is a function of camber as well 
as thickness, but since airfoils suitable for use in lift- 
ing airscrews must have small cambers to prevent undue 
blade twist and possible flutter, we may neglect this ef- 
fect. 

Fig ur e -of -merit curyegt^ A study of the effects of 
varying the different parameters can bo made with the aid 
of curves of figure of merit M, plotted against blade- 
incidence factor 6^. Recalling that 
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T 3/2 



2 Q CT 

, m 3 /2 

_ 1 ZSL 

4a 

we note that M = f (e, -^j, Bg.) , assuming variations in 

a 0 to "be negligible. 

The effect of € upon M for the constant-incidence 
rotor may he judged from figure 16, which was computed from 
the estimated full-scale data of table VII. Values of t 
of 0.1, 0.3, and 0.5 and a range of G a from zero to 9 are 
shown, and median values of the other parameters were used. 
In this and subsequent figures the curve of ideal figure of 
merit Mj , is included for comparison. 

It is interesting to note that e has practically no 
effect up to % = 1.5. However, at 6j = 4.8, which is 

the upper limit of a 0 = 5.75, for a = 0.06 an apprecia- 
ble divergence is evident. The curve ( = 0.1 will doubt- 
less never be reached unless very much larger rotors than 
now exist are developed, whereas e = 0.5 is representa- 
tive of small rotors such as might be used as test models 
in a large wind tunnel. 

In figure 17 similar curves are shown for 8 = 0.006, 
0.009, and 0.012, computed for the median full-scale value 
e = 0.3. Here the lowest value of . 6 represents a large 
rotor with a thin blade profile, and the highest, a small 
rotor with a thick profile. The value 6 = 0.009 repre- 
sents the probable full-scale median. Comparison of fig- 
ures 16 and 17 shows that for 6 a > 3.5, € has a greater 
effect on M than 5, whereas for 6 a < 3.5, the reverse 
is t rue . 

An idea of the effect of blade twist on M can be ob- 
tained from figure 18 for the two cases 6 = constant and 
6 

6 - -~ (constant pitch). Since no suitable experimental 

data are available for the constant-pitch airscrew of con- 
stant chord and profile, it was necessary to compute the 
curves of figure 18 from the theoretical thrust and torque 
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equations. The relativ e values of M for the two types 
of "blade should "be given reliably "by these curves. It 
will "be noted that the percentage increase in M of the 
constant-pitch rotor over that for the constant-incidence 
rotor is greater for small values of 6 a and the ratio of 

the values of M max is about 1.07. 

The theoretical ideal figure of merit Mj , for the 

constant-pitch rotor has already "been shown to he a con- 
stant and equal to unity. For the constant-incidence ro- 
tor, M j , is not exactly constant, "but increases slowly 
with Gqt . This variation is so small, however, that Mj = 
0.94 = constant may he considered a satisfactory repre- 
sentation. The curve of Mj calculated from the modified 
data of tahlc VII is included in the figure and shows. a 
large deviation from this theoretical value at small val- 
ues of 6 a . 

In connection with figure 18 it should, of course, 
he kept in mind that the curve p = constant represents a 
different blade twist for each value of 6 a . If the "blades 

were merely turn.ed as a whole without change of twist, the 
values of M would lie somewhcire "between the two curves. 



The foregoing analysis of the effect of changing the 
"blade parameters clarifies somewhat this phase of the lift- 
ing airscrew problem. However, there still remains a need, 
for more consistent and extensive large-scale data on the 
characteristics of airfoils suitable for lifting airscrews. 

Having established full-scale values of 'the parameters, 
it is now possible to construct a general f igure-of-merit 
qhart from which, together with figure 15, the hovering 
characteristics of any constant-chord and incidence-lifting 
airscrew can be quickly derived. This chart is shown in 
figure 19 and is based on the following parametric values: 



5 . 75 



£ 



0.3 



8 



0.006 and 0.012 



6 



0 



to 



3.2 
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The narrow cross-hatched Land in this figure represents 
the region of l/L««, (limit of a n = 5.75) for different val 

ues of 8 and a. This region indicates the significant 
fact that the "lifting efficiency" of all rotors of this 
type will "be very nearly 0.80 with extreme deviations of 
±0 # 04. The most efficient rotor for hovering flight should 
ohviously have a large solidity and a small -"blade profile , 
drag, "but the small solidity rotors essential for efficient 
forv/ard motion would not "be greatly inferior when hovering. 

As a matter of theoretical interest the curves of Mj 

and = 0 are included in this chart, and all the curves 

a G 3 

are carried up to 6g = 9 for the same reason. In this 

connection it will "be noted that the full-scale limit of 
a Q = 5,75, i.e., a = 14°, enahles the maximum value of 

M to "be very nearly reached in most cases. 



ZOVERING-FLIGHT PERFORMANCE 



The ho vering-f light performance characteristics of any 
helicopter may "be readily determined "by the simple method 
developed hy Glauert in reference 2. This method has "been 
somewhat modified by the writers and has "been reduced to 
two simple charts from which altitude performance can "be 
easily calculated. 

The figure of merit may "be expressed as 



1 1 / 2T T a 



3 /2 



M = — — — /_ —g- . Jl (19) 

2 OQ J p tt R %q 

where OQ is the power required. And "by analogy we may 
write : 

i B ii rriTZ (20) 

2 P I J p TT R 

whore P 1 = T| ? is the power available at rotor 
P, power of motor 
T| f efficiency of drive system 



and 



W, weight 
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The factor T) is very nearly unity for ordinary pro- 
pulsive airscrews, but for the lifting airscrew it will 
vary considerably, depending upon the type of drive used. 

From the above expressions, it is at once apparent 
that for hovering flight at a given altitude 

| = II 

5 QQ, 



which is the ratio of power available to power required, 
and that hovering flight is therefore Dossible only when 
u ?». 

The expression for M also shows that for a given 
thrust, the power required will be inversely proportional 
to the radius R of the rotor and, consequently, for max- 
imum economy the rotor should be as large as possible. 

Using the subscript zero to denote sea-level condi- 
tions, 

* = l-f / P -0 = f (h) f (H ) (21) 



and adopting the relationships for -oower and density varia- 
tion with altitude given by Diehl in reference 8, figure 
20 has been prepared by plotting M/N 0 vs. H. 

Since for hovering flight M = K, the following rela- 
tionships at sea level may be obtained from equation (20): 




In this equation two loading factors are apparent and these 
will be designated as 



| ) ss l the "ideal effective" power loading 

\MP 1 Q > P 



( W 

( — g- J a 1 , the rotor disk loading 
v v,H ' 



Thus equation (22) becomes 
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(23) 



and this relationship is plotted in figure 21 in terms of 
e n gin eering units. 

The concept of an "ideal effective' 1 power loading is 
useful in that it makes i^ossible the simple relationship 
of equation (23) and the single curve of figure 21, which 
is perfectly general. The form of this fictitious power 
loading shows that it merely signifies a higher value for 
actual rotors and, therefore, an inferior performance as 
compared with the ideal case of Mj = 1, 

The power ratio for hovering flight may now "be ex- 
pressed in terms of the loading factors as follows: 



JL 

h J 1 s 



(24) 



Determination of the rotor tip speed can he made quite 
simply. If we write that 

3 /2 



_T_ 
P ' 



and 



'MP 

out remembering that M 



JL _ a V 
pt ~ n r q a 



a T 



a 



M Q R Q 



1 - 



m 3 /2 

, we obtain 



2 Q CT 



Q R 



a 1 B T a 



(25) 



and since for a given rotor and blade angle, 



and To- 



are constant, the tip speed is inversely proportional to the 
"ie'eal effective" power loading, 

^i-iii^d-e^^G rformari.ee .._calculat_i on^, - The following sim- 
ple expressions, converted to engineering units, enable the 
calculation of the performance parameters by means of which, 
with the curves of figures 20 and 2"1 , the hovering perform- 
ance characteristics of any helicopter may be determined. 
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Lp = [.•■"■) in lb. ner hn . 

L „ = (— — -g ^ in It), ner sq.ft. 



(A) 
(3) 



M 38 (r \ 

1 o L p L s 

QR = in ft. ver sec. (D) 

LpCcr T ff l/S 

The following example is included to show how the 
method is used. Since adequate data on recent constant- 
chord helicopters are not available, the C 30 autogiro has 
"been selected for this purpose. This machine has the fol- 
lowing approximate characteristics : 

Rotor radius, R 18.5 ft. 

Normal loaded weight, W 1,800 lb. 

Assumed effective power, T 1 . . . . 120 hp. 

Solidity, p 0.05 

Minimum profile-drag coefficient, 8 0.008 

1. The value of M can he determined from figure 19, 

since — = 0,8, and is seen to "be 0.81 at a blade inci- 
4 a 2 

dence factor of 6^ = 7 or a hlade incidence of 

G =s 57.3 a e a 
= 20.1° 

2. The sea-level minimum loading factors will "be 

/ 1800 \ 

Lp s 1 = 18.5 Id . ner hp. 

V M3.81 x 120^ 

iSOO 

L g = g" = 1 • o 8 lb. per so. it. 

l - TT 18.5 
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3. The power ratio ' from equation (C) is 
M 38 • 



n o 18.5 x (1.68) 



1.59 



4. The absolute ceiling for hovering flight can he 
obtained immediately from figure 20, as indicated by the 
dotted lines, and is 

H r a 8,700 feet for — ■ 1.59 

N o 

5, The maximum load that can "be sustained at sea level 
(neglecting ground effect) is obtained by noting in figure 
21 that the upper limit of disk loading L M , which occurs 

at L«p s 13.5, is 4.2 and, therefore, 



- l 

L 

1.68 



s -i 



6. 

1 ecting 
s imilar 
1 oa&ing 



The minimum rotor pow.er required at sea level (neg- 

also obtained from figure 21 in a 
that the upper limit of power 

L a s 1«6S, is 29.3, and 



ground effect) is 
fashion by noting 



which occurs at 



therefore 



18 5 

Kin = 7~ ? nax = 120 = 76 ^ • 

J p C_ . o 

x o 



7. The rotor tip speeds corresponding to the maximum 
and minimum loadings are obtained from equation (D). The 
thrust coefficient must first be obtained from figure 15 
and for a = 0.05 Is found to be 9.0. Thus, for 

18.5 corresponding to ^ max = 4,500 pounds. 

1 100 

Q.R = — — m = Z?6 ft. TDer sec. 

18.5 X 0.05 X 9f 

and for L_ Q s 29.3 corresponding to ^ n ± n 55 1 , 800 pounds. 
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CONCLUSIONS 



The following conclusions can "be drawn from the re- 
sults of this investigation on the lifting airscrew: 

1. The vortex concept used in deriving the expression 
for the axial induced velocity yields results identical 
with those obtained "by the use of the momentum equation. 

2 # The experiments amply verify the theoretical con- 
clusion that the solidity can he eliminated as a separate 
parameter insofar as thrust and torque are concerned if the 
minimum torque is subtracted from the total torque. 

3. The experiments verify the assumption that for a 
given solidity the results are independent of the number of 
blade s • 

4. The theoretical equations disclose the existence 
of certain airfoil parameters with the aid of which scale 
effect can be accounted for quite simply. 

5« General full-scale curves of thrust and torque for 
the constant chord and incidence rotor can be estimated 
with reasonable confidence by a judicious combination of 
the theoretical and experimental data of this report. 

6. In general, it appears that the optimum figure of 
merit for the full-scale constant chord and incidence rotor 
may be taken as 0.80, and for the constant chord and pitch 
rotor as 0.85 with a probable maximum error for extreme 
cases of not more than i0.04. 

7. Helicopter ho vering-f light performance can be rapid- 
ly determined by means of the two simple ceiling and loading 
charts of this report when the figure of merit and thrust 
coefficient are known. 
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TABLE I 

Helicopter Model Tests - Static Thrust 
2-"blade rotor, <J = 0.0424 
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TABLE II 

Helicopter Model Tests - Static Thrust 
3-blade rotor, a = 0,0636 
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TABLE III 

Helicopter Model Testa - Static Thrust 
4-"blade rotor, a =s 0.0849 
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TABLE IV 

Helicopter Model Tests ~ Static Thrust 
5~"blade rotorf, (J = 0.1061 
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TABLE V 

N.A.C.A. 0015 Profile Characteristics 

6- by 72-inch airfoil H.N. = 242,000 

(Position A ) 
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TABLE VI 

N.A.C.A. 0015 Profile Characteristics 
5- by 72-inch airfoil R.N. = 242,000 

(Position B) 
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Eelicopter Static Thrust 
Full-scale data (est.) 
6 = const., 3 = co, a Q = 5.75 
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Points taken from faired experimental data on four rotors» 
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TABLE VIII 

Helicopter Theory - Static Thrust 
t 

6 a const., B = oo, a 0 = 5.75 
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TABLE IX 

Helicopter Theory - Static Thrust 
p = const. , S = oo, a 0 = 5.75 
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Fig. 2 
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figure 2. 



N.A.C.A. Technical Bote Ho. 626 



Figs. 4, 5,6 




Figure 4.- Three-b laded rotor mounted Figure 5. - Rotor model mounted on blade 
on wind tunnel balance. incidence Jig. 




Figure 6.- Upstream view of airfoil model mounted for force teste in 
9 ft. wind tunnel. 
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